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The contact angle that a liquid drop makes on a soft substrate does not obey the classical Young's 
relation, since the solid is deformed elastically by the action of the capillary forces. The finite 
elasticity of the solid also renders the contact angles different from that predicted by Neumann's 
law, which applies when the drop is floating on another liquid. Here we derive an elasto-capillary 
model for contact angles on a soft solid, by coupling a mean-field model for the molecular interactions 
to elasticity. We demonstrate that the limit of vanishing elastic modulus yields Neumann's law or a 
slight variation thereof, depending on the force transmission in the solid surface layer. The change 
in contact angle from the rigid limit (Young) to the soft limit (Neumann) appears when the length 
scale defined by the ratio of surface tension to elastic modulus j/E reaches a few molecular sizes. 



The wetting of liquid drops on dcformable solids is im- 
portant in many circumstances, with examples from bi- 
ology to microfluidic devices [1-4]. When the solid is 
soft or flexible, the shape of both the solid and the liq- 
uid are determined by elasto-capillary interactions, i.e. 
by the elastic response to the capillary forces [5, 6]. Till 
date, however, the most basic characterization of wetting 
has remained elusive for highly deformable solids [7-10]: 
What is the contact angle that a liquid makes on a soft 
solid? 

The geometry of the interfaces near the three-phase 
contact line is governed by two classical laws that de- 
scribe the macroscopic boundary condition for the con- 
tact angles [11]. Young's law applies in the case where the 
substrate is perfectly rigid, with elastic modulus E = oo, 
while Neumann's law holds for liquid lenses floating on 
another liquid "substrate". A question that naturally 
arises is whether the contact angles vary from "Young" 
to "Neumann" upon reducing the elastic modulus of the 
substrate: in other words, does one recover Neumann's 
angles in the limit E — > 0? Interestingly, the ratio of 
liquid-vapor surface tension 7 to elastic modulus E has 
the dimension of a length. It has remained an object of 
discussion whether, for the solid to become highly dc- 
formable, this elastic length -f/E should be comparable 
to a molecular size [12, 13] or to a macroscopic length 
such as the size of the drop [14, 15]. 

The difficulty of the problem results from its inherently 
multi-scale nature. On one hand, the capillary forces are 
localized in the vicinity of the contact-line. On the other 
hand, the Green function giving the surface displacement 
Sh(x), induced by a Dirac force distribution of resultant 
f z applied at the boundary of a two-dimensional elastic 
medium, scales as [16] 

6h{x) ~ in\x\, (1) 
and is therefore singular at both small and large distance 



x from the contact line. An outer cut-off is naturally 
provided by the thickness h of the elastic film or the 
size of the drop [12, 17, 18]. The inner regularization 
is commonly assumed to originate from the finite range 
of intermolccular capillary forces [7, 10, 13], or by the 
breakdown of linear elasticity [19]. Hence, the transition 
from Young's to Neumann's contact angles calls for an 
unprecedented, fully self-consistent treatment of clastic 
and capillary interactions. 

In this Letter, we solve the elasto-capillary contact an- 
gle selection within the framework of the Density Func- 
tional Theory, using the sharp-kink approximation. The 
evolution of the contact angles with stiffness is summa- 
rized in Fig. 1. The central result is that the liquid con- 
tact angle is selected at the molecular scale a and there- 
fore exhibits a transition from "Young" to "Neumann" 
around a dimcnsionless number j/(Ea) of order unity. 
We propose an analytical description of this transition, 
which agrees quantitatively with the full numerical solu- 
tion of the coupled DFT and elasticity models. Above 
this transition, the elastic solid is deformed by the cap- 
illary forces over the length "f/E. When the latter be- 
comes larger than the system size (the layer thickness h 
in Fig. lc), the elastic deformation saturates. 

Density Functional Theory - The multi-scale nature of 
elasticity makes it convenient to treat the wetting inter- 
actions in a continuum framework, such as the DFT. We 
consider a simplified DFT model in which the solid and 
the liquid are treated as homogeneous phases that mutu- 
ally attract, while the interface is assumed to be infinitely 
thin [20-22]. This model captures the microscopic prop- 
erties such as the stress-anisotropy near the interface, the 
disjoining pressure and the line tension and is consistent 
with macroscopic thermodynamics in the form of Laplace 
pressure and Young's law [20-23]. 

The idea underlying this DFT model is to separate 
the molecular interactions into a long-range attractive 
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FIG. 1: Geometry near the three-phase contact line obtained by coupling elasticity to a DFT model. Contact angles continuously 
vary from Young's law to Neumann's law by reducing the stiffness of the solid, (a) Rigid solid, j/(Ea) <C 1. The surface is 
undcformed and the liquid contact angle follows Young's law down to molecular scale a. (b) Soft solid, ^/(Ea) 3> 1. Surface 
elasticity is negligible on the scale of molecular interactions, and the contact angles obey Neumann's law. The solid is deformed 
over a distance ~ j/E from the contact line, (c) Very soft solid, "//(Eh) ~ 1. The change of the contact angles saturates when 
y/E becomes comparable to the thickness of the elastic film. The solid angle measured at scale h becomes identical to the 
microscopic solid angle at scale a. 



potential f(r), which takes into account the pair corre- 
lation function, and a hard core repulsion that acts as 
a contact force. For van der Waals interactions this po- 
tential decays as 1/r 6 , which is cut-off at a microscopic 
distance r — a that corresponds to the repulsive core. In 
the model, it turns out that all the capillary forces can 
be expressed in terms of the integrated potential [13], 

$ a/3 (r) = p a p p / dr' (p a p(\r-r'\). (2) 

J a 

This represents the potential energy in phase /3 due 
to phase a, where the phases can be liquid (L), solid 
(S) or vapor (V). p a and pp are the corresponding ho- 
mogeneous densities. The repulsive core at r = a en- 
sures that the integrals over the entire domain a con- 
verge, and is modeled by an isotropic internal pressure 
that ensures incompressibility. As detailed in [13, 24], 
the model distinguishes three types of attractive inter- 
actions: liquid-liquid, solid-solid and solid-liquid inter- 
actions, which can be expressed directly in terms of the 
surface tensions 7, 75^ and jsl [13, 21, 25]. The liquid- 
vapor surface tension 7 characterises the liquid-liquid in- 
teractions. The strength of the solid-liquid interactions 
is characterized by Young's contact angle 9y, defined by 
cos 9y — ("fsv — 1sl)/i- The interaction with vapor can 
be neglected in the limit of a low vapor density. In the 
full DFT numerical calculation, the equilibrium shape of 
the liquid-vapor interface is determined iteratively using 
the procedure described in previous papers [22, 23]. 

Selection of the liquid angle - An important feature is 
that the strength of the capillary interactions depends 
on the geometry of the deformable solid. We consider 
the reference case of a solid shaped like a wedge of angle 
9s (upper inset of Fig. 2). Similar to the case of a fiat 
surface, the force acting on a corner of liquid depends 
only on its angle 6l at a large distance from the contact 
line, and can be determined exactly by integrating over 



all the interactions in the DFT model [21, 22]. This force 
on the liquid corner consists of three contributions that 
are sketched in the lower inset of Fig. 2: (i) the force ex- 
erted by the solid (solid- liquid interactions, black arrow), 
(ii) the attractive force exerted by the rest of the liquid 
(liquid- liquid interactions, white arrows), and (iii) the re- 
pulsive force exerted by the rest of the liquid, induced by 
the presence of the solid [26] (liquid- liquid interactions, 
red arrow). This last force arises because the presence 
of the solid leads to an increase of the liquid internal 
pressure near the solid-liquid interface. 

The balance of forces in Fig. 2 provides the equilib- 
rium Ql for arbitrary 9s (details are worked out in the 
Supplementary [24]): 

cos 6l = ^ [cos 9y [1 — cos 9s] 

- sin 9 S \ - ' 2 — r , cos 2 6>y . (3) 

V 1 - cos 9 s 

This result is independent of the microscopic length a 
and the functional form of ip(r). For a flat surface 
(9s = 7r), the solid-on-liquid force is oriented vertically, 
with fsL = 7sin#L. In this case the force balance re- 
duces to Young's law, and the liquid angle 9l = Oy. 
However, (3) predicts that 9^ increases when #5 is re- 
duced (Fig. 2, solid line). Physically, this is due to the 
reduction of the solid volume for smaller #5: this lowers 
the total solid-liquid interaction, making the solid wedge 
more "hydrophobic" . 

Selection of the solid angle - If the phase S behaves as 
a perfect liquid, its mechanical equilibrium gives a sec- 
ond equation for the angles. This can be deduced from 
(3) by exchanging the roles of L and S, which indeed 
result into 9s and 9l according to Neumann's law [24]. 
In the elasto-capillary problem, by contrast, the solid S 
can resist shear. One therefore needs to express how the 



3vr/4- 



tt/2- 



tt/4- 





/ 


Os 






tt/4 



tt/2 



3tt/4 



0S 



i 

7T 



FIG. 2: (Color online) Main graph: Relation between 8^ and 
8s predicted by the DFT model. The solid line is the analytic 
formula (3) for 8y = 0.96. Symbols are the angles obtained 
numerically for the normal force transmission model (■) and 
the vectorial force transmission model (•), as defined in the 
text. Upper inset: definition of 8l and 8$- Lower inset: forces 
acting on a corner of liquid (bright region, light blue). Black: 
force exerted by the solid. Red: repulsive liquid-liquid force 
induced by the presence of the solid. White: attractive force 
exerted by the liquid, due to the missing half domain of liquid. 

capillary stress cr applied at the free surface deforms the 
solid. We treat the substrate as an incompressible elastic 
body (Poisson ratio v = 1/2) with Young's modulus E, 
as is typical for soft elastomers. Introducing the Green's 
function R, which depends on the elastic properties and 
the geometry of the substrate, we get the surface dis- 
placement: 



Sh(x) 



1 f°° 

— / H(x — x ; h) ■ cr(x )dx . 

& J-oo 



The contact line is considered to be invariant in one 
direction, so that R and cr have two components cor- 
responding to the normal and the tangential directions 
to the substrate. The elastic kernel requires a cut-off 
length at large scale, which for our numerical calcula- 
tions arises due to the finite elastic film thickness h [14]. 
The capillarity-induced cr can be expressed in terms of 
the $ Q) g [24], and the integrals of (2) can be evaluated 
numerically for arbitrary shape of the liquid and solid do- 
mains. This closes the elasto-capillary problem and the 
resulting numerical profiles are provided in Fig. 1. 

At intermediate distances from the contact line, a <C 
x -C h, the Green's function for the elastic response is 
given by Eq. (1). The slope of the solid-liquid interface 
thus scales as Shf ~ f z /(Ex). Importantly, the angle 8l 
of the liquid is selected at the micro-scale a. Therefore, 




FIG. 3: (Color online) Forces acting on the corner of solid near 
the contact line (indicated by the bright (light orange) region 
near the contact line), (a) Normal force transmission model. 
Black: force exerted by the liquid. Red: force exerted by 
the solid due to pressure-build induced by the liquid. White: 
force exerted by the solid, due to the missing half domain of 
solid, (b) Vectorial force transmission model. The difference 
with respect to (a) is the absence of pressure build-up in the 
surface layer (red). 

the relevant solid angle 8s induced by elastic deforma- 
tions must be denned at that scale. This is confirmed 
by the agreement between the prediction of (3) and the 
numerical solution of the fully coupled elasticity-DFT 
model: the symbols in Fig. 2 are obtained by measur- 
ing 9s in the numerics at a distance a from the contact 
line. With this information, one can obtain an approx- 
imate equation for the selection of 8s by evaluating (1) 
at x = a: 



5h' ~ tan 



7T - V S 



A 

Ea 



(4) 



The force acting on the solid corner. - The final step 
is to express the vertical force f z exerted on the solid 
corner in the vicinity of the contact line (bright, light 
orange region in Fig. 3). Using the approximation that 
the solid domain is a perfect wedge and assuming that the 
liquid is at equilibrium, we can derive the tangential and 
normal components of this force due to the liquid-solid 
interaction [241, 
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fl-COSfly) + 8 L 
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(5) 



(6) 



As emphasized in recent papers, this force is oriented 
towards the interior of the liquid and therefore presents 
a large tangential component, even in the limit of a flat 
surface [13, 27]. 

To express the solid-solid interactions, we need to 
model the mechanical behavior of the surface layer of 
the substrate. We consider two extreme cases of how the 
liquid-on-solid force can be transmitted to the bulk of the 
elastic solid. First, one can assume that only the normal 
stress is transmitted, as would be the case for a liquid. 
In terms of forces on the solid corner (bright region in 
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FIG. 4: Transition of the contact angles 9l (white), 8s (gray) 
upon increasing the "softness" parameter 'yf(Ea) for 9y = 
0.96 and *ysv ~ "/■ Symbols correspond to DFT numerical 
solutions for h/a — 1000 for the normal transmission force 
(□) and for the vectorial transmission force (o). 

Fig. 3), the tangential component of Jls is balanced by 
a pressure build-up in the surface layer. This is repre- 
sented by the red arrow in Fig. 3a (in perfect analogy 
to the red arrow in the liquid in Fig. 2). In this case of 
normal force transmission, the total vertical force reads 

f» = f2s cos ^SL+(-fis+l-lSL) sin6>s L -7 S v sin6»sy . 

(7) 

The angles 9sl and 9sv are defined with respect to the 
undisturbed solid surface (Fig. 2). 

Alternatively, one can hypothesize a perfect vectorial 
force transmission, for which there is no such pressure 
build-up in the surface layer (Fig. 3b). We recently pro- 
posed an experimental test aiming to discriminate be- 
tween the two force transmission models: it turned out 
that the vectorial transmission model is the correct de- 
scription for an elastomer [27] . Then, the tangential force 
exerted by the liquid is transmitted to the bulk of the 
elastic body, and the total force on the solid corner be- 
comes (Fig. 3b) 

fz = fhs cos0sl + (-/is - isv) s\tl9 S l - 7sv sm0 S y . 

(8) 

Discussion - The system of equations (3-6), closed by 
Eq. (7) or Eq. (8), gives a prediction for 9 L and 9s, and 
forms the central result of our paper. It involves three 
dimensionless parameters: the elasto-capillary number 
j/(Ea) and the surface tension ratios 7sy/7, 7sl/7- The 
resulting contact angles are shown as solid lines in Fig. 4. 
The symbols provide a direct comparison with numerical 



solution of the full elasticity-DFT model. 

In the limit of strong elasticity, -//(Ea) <C 1, one re- 
covers Young's angle 9l = 9 Y and an undeformed solid 
9$ = it, regardless of the model of force transmission. In 
the opposite limit of a soft solid, j/(Ea) 1, the elas- 
ticity is too weak to resist any force near the contact line. 
In this limit one thus finds that f z = 0, corresponding 
to a perfect balance of capillary forces. Elementary alge- 
bra shows that for the normal force transmission model 
(Eq. 7), the angles 9l,9s perfectly correspond to Neu- 
mann's triangle in the soft limit. Treating the capillary 
force as perfectly localized [15], i.e. a — 0, corresponds 
to this limiting case. However, a continuous transition 
between "Young" and "Neumann" is observed in Fig. 4, 
with a crossover around j/(Ea) ~ 1. A similar transition 
of the contact angles is observed for the vectorial force 
transmission model (Eq. 8), but with angles in the soft 
limit different from Neumann's law. 

Capillarity induced surface deformations are known to 
have major consequences for condensation of drops on 
soft solids [28] , as used for micropatterning of polymeric 
surfaces [29] , mechanical stability of gels [30] , or wetting 
on very soft coatings. The presented generalization of 
the classical laws for the contact angle will provide the 
essential boundary condition for such problems. 
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In this Supplementary Online Material, we describe the technical details ol the results provided 
in the main paper. First, we discuss how to compute capillary forces in the vicinity of the contact 
line within the Density Functional Theory, in case the solid and the liquid phases can be considered 
as perfect wedges. In this framework, we determine the equilibrium of a liquid when the solid is 
not a flat substrate, but a wedge. We show how Young's law and Neumann's law are recovered in 
particular limits. Second, we present a macroscopic linear model for the elastic response of the solid 
to capillary forces. This is the basis for the analytical and numerical results presented in the main 
paper. 



CALCULATION OF CAPILLARY FORCES IN WEDGES 



The DFT framework in which we perform our analysis separates the molecular interactions into a long-range 
attractive potential and a short-range hard core repulsion [1, 2]. The short-range repulsion is described by an internal 
repulsive pressure p r that ensures incompressibility of both the liquid and the solid - this is a reasonable approximation 
for elastomers. This pressure acts as a purely normal contact force, and is continuous across the interfaces. Using 
Eq. 2 of the main text to define 3>sl and the mechanical equilibrium allows one to express p r in the liquid as 

(see [2] for the derivation): 

Pr + + <&ll = est . (1) 

The force Fls exerted by a wedge of liquid on a subsystem of the solid of volume Vs, can be expressed as the 
resultant of the volumetric attraction and the repulsion at the liquid-solid interface Sls- Within the same framework, 
one can compute the force F$l exerted by a wedge of solid on a subsystem of the liquid of volume Vl- One writes: 

Fls=-[ V$ LS (r)d 3 r- f p r nd 2 r , (2) 



F S l=- V$sL(r)d d r- / p r nd 2 r . (3) 

JVl JSls 

When the liquid and the solid domains are perfect wedges, these integrals can be evaluated explicitly, and expressed 
in terms of the surface tensions 7, Jsl, isv- This will be derived below and lead to the equations of the main paper, 
in particular Eqs. 3, 5 and 6. 
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FIG. 1: Definition of the attractive potential $^2 (z, a), created by a wedge (of angle a) of phase 1 on a point of phase 2. The 
expression is given in (4). 



First, we determine the equilibrium shape of a liquid wedge on a given solid wedge, and the force exerted by 
the liquid on the solid in the vicinity of the contact line. This calculation can be reduced entirely in terms of the 
interactions between wedges, as in Fig. 1. The attractive potential $]^(z, a) created by a wedge of a phase 1, and 
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characterized by an angle < a < tt /2, on a point of phase 2, at a distance z from the corner, and aligned with one 
of the two edges (see Fig. 1) can be analytically expressed as: 



(z, a) = J r a — arcsin ^ - sin <pi2(r)dr . (4) 



where (p\2{r) = p\Pi 1°°^ { P\i{\/t' 1 + y 2 )dy is the potential created by an infinite line at a distance r from it. Note 
that by construction, this function must be even with respect to the angle. The more general case when a point is not 
aligned with an edge of the wedge can be easily deduced from the above case by adding (or subtracting) two different 
wedges. Furthermore, as the surface tensions are calculated from the interactions between a column and half a space, 
we can obtain certain important relations governing ^12 (r) [3]: 



f 

Jo 



r (p LL (r)dr = - 7 , (5) 



°° 2 ~ / \ 7 f°° 2- , s A 7 + 7sv-7sl l + cos0y 

rip LS {r)dr = J r tp SL {r)dr = = -7 . (6) 

Here Qy is the Young' angle. These two integrals can now be used in (4) to obtain: 

f°° *w 1 \ j 1-cosa 

r*YsMdz= r^ L (z, a) dz= -? 1+ 7 9y1 27 - ^ 

J a Jo z sin a 

The attractive potential $12 (2) created by a half-space of a phase 1, on a point at a distance z from its flat boundary 
can be expressed independently or by using two wedges of 90° : 

$f 2 (z)-2^ °rarccos^ 12 (r) ( ir = 2$^(z,|) . (9) 

This quantity (noted n^z) in [2]) can be interpreted as a disjoining pressure. Furthermore, we have some useful 
symmetry properties that we can easily interpret from geometrical arguments (here we dropped the subscripts to 
clarify the notations): 

$ ff (z<0)=2$ ff (0)-$ H (-z), (10) 
<S> w (z,a > tt/2) =$ H (zsina)-$ w/ (z,7r-a) , (11) 
$ M/ (z<0,a) = $ H (0)-$ w/ (-z, 7 r-a). (12) 



Force exerted by the solid on the liquid near the contact line 



We denote by Slv, Ssv and Ssl the liquid- vapor, solid- vapor and solid- liquid interfaces, respectively. The subsys- 
tem of liquid subject to forces is bounded on the liquid side by a surface perpendicular to the liquid-solid interface, 
denoted Sll ( see Fig. 2). The angle of the liquid wedge is < 9 L < tt/2 and the angle of the solid wedge is < 9 S < tt. 
The force exerted by the solid wedge on the liquid subsystem in the vicinity of the corner is given by (3). Using Os- 
trogradsky's theorem, the volume integral in (3) reduces to a surface integral, in terms of the surfaces bounding the 
volume of the liquid subsystem: 

F SL = - f 4l (- sin e L u x + cos 9 L u z ) dS - f u x dS+ f Uf L + p ( ^) u z dS , (13) 

where u x and u z are the tangential and the normal unit vectors to the solid-liquid interface, both directed toward 
the interior of the liquid. We thus have to calculate the three different attractive potentials <J>g£ created by the 

(3) 

solid wedge on different locations of the liquid wedge, and the repulsive pressure p]~ in the liquid at the liquid-solid 
interface. Let us first concentrate on the liquid at the liquid-solid interface. The liquid is at equilibrium, so we can 
use Eq. (1) to obtain (in the entire liquid) p r + &sl + ^ll = $ffz,(0), which can be simplified to obtain: 

P? ) +* { H = *Il{x,'k-9l). (14) 
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Then, we can express the solid-liquid attractive potentials $5 2 and $3/ respectively far from the contact line at a 
distance z from the solid-liquid interface, and along the liquid-vapor interface at a distance s from the contact line: 



(2) 
SL 



4l = - 0L) + + - 7T) 



(15) 
(16) 



Using the different expressions of the potentials above (14-16) and the integrals (7,8), we can finally explicitly express 
the force (13) per unit L of contact line along the two directions normal and tangential to the solid-liquid interface: 



ft F SL ■ u x 

Jsl = £ = 7 1 1 



+ cos By) 
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= 7 



(1 + COS0y) 



sin 9l 
2 

cos 6l 



1 — cos 6l 1 + cos(6l + 1 



sin 6l 
1 + cosOl 



sm{6 L + 9 S ) 
1 + cos(6» L + 



sin 



6 L ' sm(0 L + 0s) 



1 + COS 0L 

sin 6l 



(17) 
(18) 



Liquid equilibrium 

As explained in the main text, the equilibrium condition for the liquid is obtained by balancing the solid-on-liquid 
force with the forces due to liquid-liquid interactions. The system on which the forces act is taken as the liquid corner 
depicted in Fig. 2). The liquid-liquid interactions can be separated in two contributions. First, there are attractive 
liquid on liquid forces, which act at the liquid-liquid and the solid-liquid interfaces, far from the contact line. Each 
of these two forces are equal to 7 per unit contact line, and are pulling on these two corners along the liquid-vapor 
and liquid-solid interfaces. The second is due to repulsive liquid-liquid interactions, acting near the boundary normal 
to the solid-liquid interface far from the contact line. It originates from the presence of the attractive solid. The 
resulting tangential force (per length unit of contact line) is 7 + 7sv — 7sl = 7(1 + cos 9y) and is directed toward the 
interior of the subsystem [4]. The balance of these forces, tangential and normal to the interface, read: 

7(1 + cos 0^-7(1 + cos 7^ = 0, (19) 
7S in0 L +/£ L = O. (20) 

These two equations come from the projections of the gradient of a scalar: it implies that they represent a single 
independent equation that reduces to [using for instance (17) and (19)]: 

l + cos» y = (l + cos fll )(l + ^|M). (21) 
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The unique physical solution of this equation is given by: 

1 



cos V L = - 



which is Eq. (3) in the main paper. 



2 

cos 9y [1 — cos 9s] — sin 9 s \ / 7, cos 2 9y 

V 1 - COS #5 



Young's law and Neumann's law 



(22) 



One easily shows that in the case of a flat substrate, 9s = n, the tangential force f SL vanishes due to symmetry 
reasons: the solid spans an infinite half-space in the tangential direction. We then recover 8l — 9y, [cither from (19) 
or (22)], which is Young's law. In addition we note that f§ L = 7sin0y, ensuring that we have a force balance in the 
liquid subsystem in both directions. 

By assuming the substrate is a perfect liquid, we will now show that the above formulation is equivalent to Neu- 
mann's law [5]. Denoting this liquid substrate as phase 2, and the liquid wedge on top of it as phase 1, Neumann's 
law is given by the two following equations: 

7i sin#! — 72 sin # 2 = , (23) 
7i cos 9i + 72 cos 2 + 7i2 = , (24) 

where 71 and 72 are the surface tensions of liquids 1 and 2, 712 is the interfacial tension, and (0i,02) are the positive 
angles made by the liquid wedges 1 and 2. Neumann's law is in fact ruled by two independent parameters, which we 
denote according to Young's law by: 

COS0yi = 72-712 , (25) 

71 

COS0y 2 = 71-712 . (26) 

72 

With these notations, we get 71/72 = (1 + cos 9y 2 )/(l + cos 6y\) and 712/72 = (1 — cos0yi cos 0y 2 )/(l + cos 0yi), and 
Neumann's law (23,24) can be rewritten as: 

sin 0i (1 + cos0y 2 ) - sin0 2 (l + cos0yi) = , (27) 

COS 01 (1 + COS0y 2 ) + COS 2 (1 + COS0yi) = COS 0yi COS 9y2 — 1 , (28) 

The unique physical solution can then be explicitly expressed as: 

[2 + (1 - COS0y 2 )cOS0yi] COS0yi - (1 + COS 0y 2 ) 

cos 0i = — — r . (29) 

2(1 - COS0yiCOS0y 2 ) 

or cquivalently by a function of the form 0i = N{6yi,0y2)- The second angle is given by symmetry, i.e., 02 = 
N(6 Y 2,0yi). 

Let us now show that the DFT formulation indeed reduces to Neumann's law when both the phases are liquid. For 
a given substrate wedge of angle 2 , we can see from (22) that the liquid 1 is at equilibrium when it forms a wedge of 
angle 0! = /(0y, 02)- If the phase 2 is a liquid as well, then we can use the same equation (22) with exchanged roles 
of phase 1 and 2. This gives two equations for the contact angles: 

01 =/(0ri,0 2 ) , (30) 

02 = f(0Y2,9 1 ), (31) 

where 0y2 is defined in (26). Therefore, 0i obeys the following implicit equation: 

0i= f{9 Y1 J{9 Y2l 8 1 )) . (32) 

We can easily verify that this equation has a unique solution in the interval [0,7r], which is given explicitly by (29). 
By analogy, we also recover the second angle 02 with the symmetry 9y\\9yi- 
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Force exerted by the liquid on the solid near the contact line 



Contrarily to the liquid, the capillary forces on the solid need not be at equilibrium for the solid: a resultant force 
can be balance by the elasticity. By analogy with the determination of the solid on liquid force, we compute the 
capillary force (2) exerted by the entire wedge of liquid to a subsystem of solid bounded by a surface perpendicular 
to the liquid-solid interface, noted Sss : 



= I *ls (- sin su x - cos 6 s u z ) dS- ( u x dS- ( Uf s + pfA 

JSsv J Sss J$sl 



u z dS . 



(33) 



The attractive potentials = $^ s (s, tt - 6 S ) + ^%{s, 9 L +0 S - n) and $^ = <§>% s {z) are analogous to what we 
obtained in (16) and (15). The only fundamental difference is that for this case the pressure is determine by the liquid 
(and not the solid) equilibrium (14), so that: 



p (3) + $(3) = $W l{x ^ _ q l) + < S >W l{x ^_ 9s) _ 3>W s{x ^_ q l) 



(34) 



Therefore, after some simplifications, we can determine [from (33)] the normal and the tangential components of the 
liquid on solid force (per unit length of the contact line): 



fls — 

rn 

JLS — 



Fls ■ u x 
L 

Fls ■ u z 



1 + cos 9y 



1 + COS By 



1 — cos 8 s + sin 9 s 
sin 9 s + cos $s 



: z ~ o I. - costf.s +smtf.v — ia ~ a , ' ■ (35) 

^ - cos 9 Y 1 + cos 9 L 1 
+ — sin^J ' (36) 

These expressions can be rewritten in a more compact form as Eqs. (5) and (6) of the main paper. When #5 = tt, then 
9l = 9y [using (22)] and /£ s = 7 sin 9y [using (36)], as we would expect from the vertical action- reaction interaction 
of two full wedges, and f LS = 7(1 + cos 9y) [using (35)], as was shown in [2, 6]. 



1 + cos(fl L + 9 S ) 

sin(0 L + 9 S ) 
+ cos(6> L + 9 S ) 
sin(9 L + 9 S ) 



MACROSCOPIC ELASTO-CAPILLARY MODEL 



The description of the solid is more subtle than that of the liquid. The bulk of the elastic substrate is described as 
a homogeneous incompressible elastic medium. However, one needs to define the effective boundary conditions to be 
applied to the bulk elasticity equations to compute the deformations. As discussed in the main text, one can consider 
two models that will be referred to as normal force transmission and vectorial force transmission respectively. Below 
we present the expressions for the surface stress a s used for our numerical results. 



Curved interfaces interactions 



For later reference, we first determine the attractive potential at an interface between two phases 1 and 2 which 
presents a curvature n. For na <C 1, this potential (Eq. 2 in the main text) reduces to: 

/>oo 

*i2 = *S(0)-«/ r 2 ^ 12 (r)dr, (37) 
Jo 

where, k is defined as positive when phase 1 is smaller than phase 2. In (37), the integral term is always negative, and 
depending on whether one considers liquid-liquid, solid-solid or liquid-solid interaction, becomes equal to —7, — 7sv 
or —7(1 + cos9y)/2 [see (7,8)]. The Laplace pressure originates from this potential. 

Normal force transmission model 

An elastomer can be described as an incompressible soft solid within which the elastic stress can be written as: 

E 

a ij = -PrSij + IT (9i u J + djUi) , (38) 



G 




where Ui is the displacement vector, and E the Young's modulus. This solid is assumed to be submitted to the 
long-range attractive potential. Then it can be rigorously shown (see the proof below) that the the external solution 
must have an effective boundary condition: 

S s = ($i <? + $SS - ®SL ~ $LL + ^LL (0) - $f S (0)) ft . (39) 

Note that this is the opposite of the surface force. When the liquid phase is expressed by a vapor phase, this expression 
reduces to: 

s s = ($ss - *f s (o)) n . (40) 

The numerical solution presented in the main paper are based on numerical evaluations of the surface stress (39,40). 

The surface stress <j s acts as a generalized capillary pressure. At distances from the contact line much greater than 
the molecular cut-off a, the curvature is small so that [using (37,40)] the capillary pressure on the dry part reduces 
to the standard Laplace pressure: 

<?s = isv k n . (41) 
On the other hand, the Laplace pressure on the wet part is 

<? s = Isl K n . (42) 
Proof — The incompressibility condition reads diui = or in vector notations: 



V-u = 0. 
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We first assume that the constitutive equation holds in the interfacial zone where the elastomer is subjected to two 
long range volumetric interactions which derive from the potentials &ls an d $ss that are associated to liquid-solid 
and solid-solid interactions respectively. The equilibrium condition reads: 

-V( Pr + *is + *ss) + £V 2 S=0. (43) 

$L5 and &ss result from an integration over the volume of the liquid and the solid. They vary over a scale of few 
molecular sizes along the normal to the solid interface and attain their bulk values at distances away from the interface. 
Therefore, and $55 are only important in a thin boundary layer across the interface. We posep = p r + $LS + &SS 
and &ij = (Jij — (&ls + ®ss)dij- Then the equation reduces to: 

- Vp + EV 2 u = (44) 

in the whole domain. This is the standard elasticity problem. However the effective boundary conditions are modified. 
More precisely, there is no excess quantity in the strain tensor, but the effective stress applied to the bulk is modified. 
As the real stress cry is continuous, the pseudo-stress try is not. 



Vectorial force transmission model 



We can alternatively propose a model of elastomeric surface that allows for a fully vectorial transmission of the 
liquid-on-solid force. The polymeric chains are attracted in the direction of the liquid and thus transfer the total 
force exerted by the liquid on the solid. There is no force gradient building up in the surface layer. Integrating the 
equilibrium equation from the surface to the bulk, we see that the effective stress perfectly balances the excess forces 
due to the long range interactions: 



c?, = t ■ V 



$ LS dh i 



(45) 



+ ($LS + $ss ~ - $LL + $£ L (0) - #f s (0)) n 

This is the expression for <r s used in the numerical solution of the vectorial force transmission model. 
Again, looking at a distance from the contact line large compared to a and using 

t ■ Vt = —nn 1 

we get in the dry region a surface stress: 

<? s = (isv + 7) « n 
Beside, the Laplace pressure in the dry part is still 

<? s = isv k n. 



(46) 



(47) 
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